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Abstract 

The Jarzynski Equality relates the free energy difference between 
two equilibrium states of a system to the average of the work over all 
irreversible paths to go from one state to the other. We claim that 
the derivation of this equality is flawed, introducing an ad hoc and un- 
justified weighting factor which handles improperly the heat exchange 
with a heat bath. Therefore Liphardt et al's experiment cannot be 
viewed as a confirmation of this equality, although the numerical de- 
viations between the two are small. However, the Jarzynski Equality 
may well be a useful approximation, e.g. in measurements on single 
molecules in solution. 



1 Introduction 

Ever since Jarzynski derived a remarkable equality, the Jarzynski Equality^^'^' , 
(JE), his results have been widely used: duplicated theoretically and tested 
experimentally. It is impossible to do justice to all publications that have 
appeared in connection with the JE: a survey article would be needed for 



that and, in fact, one covering part of the material has been written by 
Ritortt^l However, in this paper we want to scrutinize the original general 
derivation of the JE^^'^J and discuss briefly a crucial experiment by Liphardt 
et al'^] which was designed to check the JE. The communities accepting the 
JE consist overwhelmingly of chemists and biophysicists to which also one of 
us (DM) belongs, while the physicists (to which EGDC belongs) have divided 
opinions. We hope that this paper may clarify and unify the various aspects 
of the JE. 

The main point of this paper is to argue that the JE is not an equality in 
any mathematical sense, but can be a useful approximate equality in certain 
important fields, like, e.g. the study of single molecules in solution. 

The JE equates the difference between the free energy of two equilib- 
rium states A and B oi a system in contact with a heat bath, with an 
average over the irreversible work done over all paths from A to B. This 
work is done in general by a set of external force centers, which can be 
characterized by a set of time dependent parameters {Aj(t)} in the Hamil- 
tonian of the system: H = H{{pi{t)}, {gi(t)}; {Aj(t)}) = H{T{t); \{t)). Here 
r(t) = {pi{t)}, {qi(t)}, defines the microscopic state of the system in its phase 
space, by giving the values of the momenta {pi{t)} and coordinates {qi{t)} 
respectively, of the i = 1,2, ...,N particles of the system at time t. For sim- 
plicity we will restrict ourselves to only one parameter A(t), so that J = 1. 
One is interested in the Helmholtz free energy difference AF = Fb — Fa 
between the initial and final equilibrium states A and B. By considering the 
work W done on the system when the external forces bring the system from 
state A to state B along all possible irreversible paths, with weights e~^^^''^. 



the JE reads: 

< e-^^"-" >= e-^^^ (1) 

where (3 = l/ksT, with T the temperature of the heat bath and ks Bohz- 
mann's constant. Since in (1) the average < > is over the work Wirr done 
over all possible irreversible paths in phase space from A to B, the deriva- 
tion of the left hand side of Eq.(l) proceeds microscopically via statistical 
mechanics. 

We will consider here the case that the system is for all times t > in 
thermal contact with a heat bath. The heat exchange with the heat bath 
can be reversible or irreversible, depending on whether it takes place with 
the internal temperature Tj of the system equal or not equal to the external 
temperature Te of the heat bath surrounding the system, respectively. In 
the former case the system will always be in a (canonical) equilibrium state 
during the work process. As is explained in textbooks of Thermodynamics, 
a slow or fast work process does not guarantee at all its reversibility or its 
irreversibility, only the equality of internal and external parameters does. 

At t = the system in state A is coupled to a heat bath of temperature 
T, so that its temperature is also T. It has then a canonical distribution 
given by'^^: 

/(r(0),0) = — e-^^(r(o)'"W) (2) 

Here Zq = Za is the canonical partition function at t = 0, where A = A(0) = 
Xa- We assume that work is done on the system during a total time r, so 
that < t < r, while A(t) goes from A(0) = A^ to A(r) = A^, the final value 
of A which corresponds to the equilibrium state B. The (mechanical) work 



W done on the system in an initial micro phase r(0) over a time t will be 

given by: 

W,.WiT.t)^r,t'?I^^^m (3) 

Jo o\ 

where A(t') = d\{t')/dt'. 

We will assume here for simplicity that the rate of change of A, i.e. A, is 
constant. 

However, if the change of A(t) with t is such that irreversible processes 
are induced in the system driving it possibly far from equilibrium, then, 
in order to reach an equilibrium final state B,X{t) must be such that the 
equilibrium state B can actually be reached. During this transition, heat 
exchange will take place (and possibly work will be done). This will depend 
on the procedure of varying A with time and on the nature of the induced 
irreversible processes. The change AF is computed by Jarzynski considering 
only the mechanical work in going from A to B. 

2 JE for thermostatted system 

Since the system is in constant contact with a heat bath, it becomes 
important whether the (constant) rate of change A of A allows thermal equi- 
librium to be maintained at all times between the system and the heat bath. 
If it does, the process is reversible^, if it does not then the process is irre- 
versible, since the internal and external variables will then not always be 



^This is usually formulated by requiring that the external changes induced in the system 
are slow when compared to the internal relaxation times needed to return to thermal 
equilibrium once a change of the system (e.g., its volume) has taken place. In that case 
the temperature of the system and that of the heat bath will always remain the same 
during a measurement. 



equal, or may not even be definable in the system. 

In addition, the reversibility can be characterized in a different way, viz., 
the strength of the coupling of the system to the heat reservoir or the rate 
of heat transfer, which we will call c. It is this quantity together with the 
work rate w = \^, which will determine whether the system will remain 
in thermal equilibrium at all times i.e., whether a reversible or irreversible 
process takes place. We can consider now several cases. 

a) The rate of work done on the system w is very small (A ~ 0) and 
the coupling c ^ also. Then depending on the ratio w/c, thermal 
equilibrium can be maintained between the system and the heat bath, 
at all times < t < r when w << c, i.e., there will be plenty of time for 
the system to exchange heat with the reservoir, so that the disturbance 
of the system due to the work done on the system and which brings 
it out of thermal equilibrium, can be readjusted for all < t < r so 
that the system is always in thermal equilibrium with the heat bath at 
temperature T. If this obtains during the entire time r that work is 
done on the system then the work process is isothermal and reversible. 

b) If the w is very small, but the coupling c is not, then the work process 
will also be isothermal and reversible. 

c) If, on the other hand, w is very large and the coupling c ^ 0, then 
there is no way to maintain thermal equilibrium during the work time 
r and the system is thermodynamically in a non-equilibrium state. 

d) However, if w is very large, then c has to be sufficiently larger to main- 



tain thermal equilibrium, if at all possible, and in general one cannot 
assume that this is so. 

3 Critique 

The above considerations are relevant because, during the evolution of an 
initial phase r(0) of the system, it will not only be subject to the mechanical 
work done by the external forces via the Hamiltonian H on the system alone, 
but also to the simultaneous energy exchange with the heat bath. Therefore, 
to know the mechanical work and the heat separately along a phase space 
path, one has, in principle, to know the microscopic state of both the system, 
the heat bath and their coupling, so as to know whether, and if so, how much 
and in what direction, heat exchange between the system and the heat bath 
has taken place. ^ 

Consequently the introduction in phase space of the weight e"'^^*'"'', for 
every irreversible (stochastic) path, where the inverse temperature (3 of the 
heat bath is used in the weight for every Wi„--, does not seem to make phys- 
ical sense. In fact, even if the mechanical work Wi„- itself could be precisely 
determined, the canonical weight e^'^'^"''' associated with it, especially when 
no internal (system) temperature is known or can even be defined - a possi- 
bility of being far from equilibrium, which Jarzynski notes himself'^^ - seems 
completely arbitrary and the use of the heat bath temperature 1//3 - the 
only known temperature available - without foundation. As a consequence. 



^For a system in contact with a heat bath the work done on the system is not just 
the mechanical work, which will change the system's internal energy (dE), but also work 
associated with the heat (energy) exchange with the bath (TdS), can already be seen in 
the isothermal reversible case, where dW = dE - dQ = dE - TdS = d{E - TS) = dF^^l 
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an average over the microscopic e"^^*'''', as carried out in Eq.(l), does not 
seem physically meaningful unless they are very close to reversible ones and 
provide then a good approximation to those. Therefore the JE is correct but 
trivial if W^ = Wrev and it seems that the use of Wirr instead of Wrev is, in 
general, unfounded. A further discussion of this point can be found in points 
of section 5. 

In the above we argued that the heat exchange has not been properly 
taken into account and an essentially only mechanical theory has been used 
to derive Eq.(l), while, however, also non- mechanical work e.g. thermal 
expansion due to heat energy has to be considered. 

A striking example of this is Jarzynski's remark'^' that it would suffice 
to let the system evolve from parameter values \a to A^* = \b regardless 
of whether B* is the equilibrium state B, whose free energy difference AF 
with state A one wants. As Jarzynski clearly states, the non- equilibrium 
state B* at Xb* = Ab, reached before the equilibrium state B, can easily 
be transformed into the desired equilibrium state B, while the amount of 
(mechanical) work done to get from Xj^ to A^ remains the same. The only 
way, however, to go from a non-equilibrium state B* with A^ to an equilib- 
rium state at B, with no more work done, is through the contact with the 
heat bath at temperature T. Since the free energy F will be a minimum in 
equilibrium, the non-equilibrium A^* = Ab state will, only if properly cho- 
sen, indeed go to the equilibrium state B, accompanied, however, with an 
exchange of an unknown amount of heat and possibly work ,which depend 
on B*, all of which will change the free energy. 

So far we have discussed the general physical theory of Jarzynski as found 
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in refs.[l,2]. This theory is supposed to hold for systems of any size. As 
Jarzynski remarkst^' to obtain observable effects the systems have to be small, 
since otherwise the fluctuations of the work values become too small to be 
observed. In this connection measurements on single molecules in aqueous 
solutions seem very appropriate as a check on the JE. 

4 Experiment 

Liphardt et al'^'^' have carried out such an experiment on the free energy 
difference (in their case the Gibbs free energy AG instead of AF) between 
the unfolded and the folded conformations of a single P5abc RNA molecule, 
suspended between two handles, in an aqueous salt solution to check the 
JE. The experiments were carried out carefully, by stretching a single RNA 
molecule many times, between the folded and unfolded conformations A and 
B. In this case the irreversibility of the procedure manifests itself in hysteresis 
curves associated with a cycle A (folded) -^ B (unfolded) -^ A (folded). A 
number of constant stretching (switching) rates of 2 - 5 pN/s (slow) and 34 
pN/s and 52 pN/s (fast), (the A above), between folded and unfolded states 
were applied and histograms were made of the work done versus the extension 
of the molecule. To obtain statistics, seven independent sets of data were 
collected for seven different RNA molecules with a slow step between two fast 
steps and about 40 unfolding-refolding cycles per molecule were performed 
i.e. about 300 independent measurements were made. Results for the average 
work differences for these three switching rates relative to the reversible work 
were plotted in bins of about 0.7 ksT . Plotting for each switching rate 



the most probable value (i.e. the maximum of the Gaussian fits to their 
histograms) of the work done on each molecule, the fast extension (switching) 
rates produced work values ^ 2 — SksT, above the expected thermodynamic 
free energy difference of about 60 ksT, because of irreversible contributions 
(apart from measurement errors) to the work, in agreement with the Second 
Law of Thermodynamics. On the other hand, an average for each switching 
rate, using the JE Eq.(l) gave values for the expected free energy difference 
within their estimated error. 

We note that the deviations from equilibrium for even the fastest switch- 
ing time, are in fact only 5% (or about SksT), which would not be unusual 
for such measurements. In addition, the thermal equilibration time is of the 
order of picoseconds, while we estimate the structural equilibration time to 
be of the order of milliseconds'^^ and the experimental switching times of the 
order of seconds (slow) to 0.1 seconds (fast). Therefore since w < c, i.e., 
essentially reversible isothermal measurements (cases a) and b) above) were 
performed. Considering the above mentioned time scales, the entire exper- 
iment could well be very close to an isothermal one for all switching rates 
used, so that the JE equality gives results very close to those of using only 
reversible paths in phase space. ^ 

A few additional remarks on this experiment follow. 
1. The Liphardt et al. data appear to fit better to the JE for large stretching 
than to the second cumulant or the "Fluctuation Dissipation Relation (FD)" : 

■^The remark by Ritortl^l that a system may be far from equilibrium even when close 
to an isothermal process because of its small size does not seem to be correct. The small 
size will allow the fluctuations to be measurable. However, they have no bearing on the 
fast change in the control parameter, needed for deviations far from equilibrium. 



< Wirr >= Wrev + /3a^/2, wheie (T^ =< W^„ > - < Wirr >^ is the variance 
of Wirr^ with Wirr = W — Wrev = Wd, the dissipative work. This, in spite of 
the fact the system is not far from equihbrium. 

We note that in the table SI of the supplementary material their < W^ > 
(scaled with j3) for various extensions and switching rates are rather close to 
a, except < Wd > in the reversible case (very slow switch) where < Wd > 
is assumed to be zero. If one computes a^/2, using the a in the second 
column of 5*1, the resulting < Wd > is too large. If, on the other hand, one 
corrects the a's in the second column for the measurement error by a simple 
subtraction of the a's at < Wd >= 0, the estimates for < Wd > are too 
small. This bracketing of < Wd > indicates that the data are compatible 
with < Wd >= a^/2. 

2. We will now derive the JE under three assumptions: 

1) The Wirr all have the weight e"'^^"''', where P is the temperature of 
the heat bath, not the system; 

2) a Gaussian assumption is made for the measured non-equilibrium dis- 
tribution functions; 

3) a FD of < Wd >= ^ has to be used. 

Assuming that the distribution functions for Wd are all Gaussian and 
using those to fit the data, the left hand side of the JE in eq.(l) can be 



^We note that their FD is derived from the JE, thereby assuming its vahdity, which 
is supposed to be checked in the experiment. This manifests itself in the appearance 
of only work contributions to a'^, while AF should also contain, in principle, cntropic 
contributions. 
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written as: 

<g-/3H^. > = I d(3Wdexp-[f3Wd-{f3Wd-f3'ay2Y/2a'] (4) 



{J d(3W,exp -ipW, - f3'a'/2f/2a']}~' 



Carrying out the Gaussian integrals leads then to: 
< e'^^-^ > = 1 or < e-/^^^-- >= e-^^-" = e"^^^ i.e. the JE. This result 
depends crucially on the validity of < Wd >= ^. 

We emphasize that two of these three assumptions, although approximately 
correct, are not justified, as explained above, since the first two are not sat- 
isfied in non equilibrium systems. 

3. We note that the often quoted claim I^'^^'^^l that the origin of the correct- 
ness of the JE is because of the over- weighting of the negative dissipative 
work balanced by the under weighting the positive dissipative work has no 
known basis. That the JE works here, in spite of the fact that it is based on 
unfounded assumptions, is because of a felicitous accident of Gaussian statis- 
tics rather than for reasons of thermodynamics and statistical mechanics. 

5 Discussion 

1. As said before, it is clearly impossible to do justice to all the publications 
that have appeared in connection with the JE. We have confined ourselves 
here to represent the contents of two articles on the general foundations of 
the JE on Statistical Mechanics by Jarzinskit^'^^. 

2. In addition we discussed an important experiment which appears to con- 
firm the JE, since the system is in aqueous solution, close to equilibrium and 



11 



the processes are (quasi) quasi-static and therefore nearly reversible, in spite 
of some hysteresis. 

3. Although in Liphardts' experiment clear hysteresis loops can be seen, they 
occur over relatively small parts of the trajectories, which mostly appear to 
be reversible. True irreversible processes, where e.g. no temperature can 
be defined in the system, have so far not been considered experimentally 
(because of the liquid surroundings used so far). Such irreversible processes 
would question severely the assignments of a canonical weights e^^^''"'', with 
f3 the temperature of the heat bath. 

We note that although to observe the JE the systems have to be small, 
there is a gap between the experiments performed on biochemical systems (in 
aqueous solutions) and those on small physical (or chemical) systems without 
an aqueous heat bath, which have not yet been performed. 

4. There are many different theoretical models, to mimic the heat bath, 
which have been used to derive the JE. They are all stochastic and, as far 
as we are aware, all use Markovian and detailed balance properties, which 
place them, in our opinion, near predominantly reversible heat exchanges. 
We confine ourselves again to two examples. 

a) This holds for Crooks' '^^ derivation of the JE for a simple Ising model. 
He discretizes there the phase space trajectory in sequences of two sub steps, 
whereby first a control parameter A does instantaneous work W on the system 
followed by an exchange of thermal energy E with a heat bath. Either 
the first forward work step is followed by a heat exchange step which lasts 
sufficiently long compared to the characteristic relaxation times of the system 
that a state of canonical equilibrium is reached (cf. Crooks' eq.(8)) in which 
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case the process is a discrete sequence of equilibrium states and therefore a 
reversible process (like in a reversible experiment) and the JE will be obtained 
as that for a reversible process. Or the time after the first forward step is not 
sufficiently long to obtain an equilibrium state, in which case an irreversible 
process takes place, which may be close to a reversible one, but no detailed 
balance obtains at any Aj.^ 

Similarly Jarzynski's procedure applied to the Langevin Equationt^', be- 
cause of the e"^^''"'' weight assumption, is restricted to systems near equilib- 
rium, with near reversible heat exchanges. 

b) Jarzynski made numerical calculations'^^ of which we will discuss only 
that on a harmonic oscillator. He considered 10^ simulations of the work done 
on a single harmonic oscillator, whose frequency is switched from ujq = 1.0 
to uJi = 2.0 over a switching time tg- The Hamiltonian of the harmonic 
oscillator is: 

where Ux switches from ujq to uji = 2uJo. If the change from ujq to uji proceeds 
infinitely slowly (tg -^ oo) and adiabatic invariance can be applied, i.e. H/u 
= constant, then Jarzynski obtains, for a canonical distribution of initial en- 
ergies a distribution function for the work W. Using this to compute the av- 
erage of the exponential work e~^^, gives j3~^ ln{uJi/uJo) = 0.693/9"^ = AF, 
since uji/ujq is just the ratio of the canonical partition functions (cf. Jarzynski 
in ref.9 eq.(59) and fig. 3 {W^ points)). One can get an idea of the average 
value of PW for fast switchings (t^ = 1), by using a similar calculation. 

^The assumed detailed balance at every step Xi is not a direct consequence of micro- 
scopic reversibility alone, and has been proved so far only for thermal equilibrium [13,14]. 
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One finds then for tlie aritlimetic average of f3W the result (— — 1) = 1. 
Adding then Wrev = 1, a value of 2 is obtained not far from the numerical 
result (ref.7, fig. 3 {W""). However, this result depends critically on the as- 
sumption of adiabatic invariance for tg = 1, while this actually only holds 
for ts ^ oo. Moreover, if for example, the energy would change as cj^, the 
Jarzynski weighting would give 2AF and the arithmetic weighting would give 
a value 4 instead of 2 for f3W. 

Jarzynski claims that the JE holds for all systems (not only harmonic 
oscillators), while we see that the JE is critically dependent on adiabatic 
invariance for finite tg, even for the harmonic oscillator. 
5. Finally we note that eq.(l) implies that a remarkable new relation would 
follow from the JE, if it were an exact equality, viz.: 

In the literature, the eq.(6), or equivalently < e~^'^'^ >= 1, is justified by 
invoking an apparently very fortuitous general cancellation due to the weight- 
ing of Wd > and W^ < with e^^^'*. It is entirely unclear to us, how this 
can be achieved for irreversible processes in a system, for all A(t) and all 
corresponding phase space paths, using in all cases the (unconnected) heat 
bath temperature (3 in the weights. 

This would imply the equality of the average over the exponentiated work 
of all irreversible paths from A to B, to the (average over all the) reversible 
path(s) from A to B. Considering the unknown nature and wide variety of 
all irreversible paths this equality does not seem physically understandable. 
If true, it would incorporate a hitherto unknown symmetry for irreversible 

14 



processes for any switching rate, i.e. a new extension of the Second Law. To 
be sure there are special cases, hke those treated in section 4, sub 2 and in 
section 5, sub 4b, where this equahty holds under certain specific assumptions 
for irreversible processes near equilibrium. 
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